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A new approach to the electroweak properties of two- 
particle composite systems is developed. The approach is 
based on the use of the instant form of relativistic Hamilto- 
nian dynamics. The main novel feature of this approach is 
the new method of construction of the matrix element of the 
electroweak current operator. The electroweak current matrix 
element satisfies the relativistic covariance conditions and in 
the case of the electromagnetic current also the conservation 
law automatically. The properties of the system as well as the 
approximations are formulated in terms of form factors. The 
approach makes it possible to formulate relativistic impulse 
approximation in such a way that the Lorentz-covariance of 
the current is ensured. In the electromagnetic case the cur- 
rent conservation law is also ensured. Our approach gives 
good results for the pion electromagnetic form factor in the 
whole range of momentum transfers available for experiments 
at present time, as well as for lepton decay constant of pion. 

I. INTRODUCTION 

The constructing of correct quantitative methods of 
calculation for composite-particle structure is an impor- 
tant line of investigations in particle physics. In nonrel- 
ativistic dynamics there exist different correct methods 
which use model or phenomenological interaction poten- 
tials. However, in the case of high energy one needs to 
develop relativistic methods. It is worth noting that now 
the experiments on accelerators, in particular at JLab 
are performed with such an accuracy that the treat- 
ment of traditionally " nonrelativistic" systems (e.g. the 
deuteron) requires to take into account relativistic ef- 
fects. Relativistic effects are important also in the treat- 
ment of composite systems of light quarks. However, the 
relativistic treatment of hadron composite systems is a 
rather complicated problem. Let us note that the use 
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of the methods of the field theory in this case encounters 
serious difficulties. For example, it is well known that the 
perturbative QCD can not be used in the case of quark 
bound states (see, e.g., 

In the present paper we will use the relativistic con- 
stituent model which describes the hadron properties 
at quark level in terms of degrees of freedom of con- 
stituent quarks. The constituent quarks are considered 
as extended objects, the internal characteristics of which 
(MSR, anomalous magnetic moments, form factors) are 
parameters of model. As a relativistic variant of the 
constituent model we choose the method of relativistic 
Hamiltonian dynamics (RHD) (see, e.g., (|-||] and refer- 
ences therein). 

The RHD method as a relativistic theory of com- 
posite systems is based on the direct realization of the 
Poincare algebra on the set of dynamical observables on 
the Hilbert space. RHD theory of particles lie between 
local field theoretic models and nonrelativistic quantum 
mechanical models. 

Contrary to field theory, RHD is dealing with finite 
number of degrees of freedom from the very beginning. 
This is certainly a kind of a model approach. The pre- 
serving of the Poincare algebra ensures the relativistic 
invariance. So, the covariance of the description in the 
frame of RHD is due to the existence of the unique uni- 
tary representation of the inhomogeneous group SL{2, C) 
on the Hilbert space of composite system states with fi- 
nite number of degrees of freedom . 

The mathematics of RHD is similar to that of nonrel- 
ativistic quantum mechanics and permits to assimilate 
the sophisticated methods of phenomenological poten- 
tials and can be generalized to describe three or more 
particles. 

The idea of this approach — RHD — is originated by 
Dirac. In js) he considered different ways of description 
of the evolution of classical relativistic systems — differ- 
ent forms of dynamics. Dirac defined three main forms of 
dynamics: point (PF), instant (IF) and light-front (FF) 
dynamics. 

RHD is based on the simultaneous action of two fun- 



1 



damental principles: relativistic invariance and Hamilto- 
nian principle — and presents the most adequate tool to 
treat the systems with finite number of degrees of free- 
dom. 

Our aim is to construct a relativistic invariant ap- 
proach to electroweak structure of two-particle compos- 
ite systems. The main problem here is the construction 
of the current operators |P-p^ . 

It seems to us that RHD is the most adequate method 
for our purpose. The use of RHD enables one to sepa- 
rate the main degrees of freedom and thus to construct 
convenient models. 

We use one of the forms of RHD, namely a version of 
the IF. 

Our approach has a number of features that distinguish 
it from other forms of dynamics and other approaches in 
the frames of IF. 

• The electroweak current matrix element satisfies 
automatically the relativistic covariance conditions 
and in the case of the electromagnetic current also 
the conservation law. 

• We propose a modified impulse approximation 
(MIA) . It is constructed in relativistically invariant 
way. This means that our MIA does not depend on 
the choice of the coordinate frame, and this con- 
trasts principally with the "frame-dependent" im- 
pulse approximation usually used in instant form 
(IF) of dynamics. |^ 

• Our approach provides with correct and natural 
nonrelativistic limit ("the correspondence princi- 
ple" is fulfilled). 

• For composite systems (including the spin 1 case) 
the approach guarantees the uniqueness of the so- 
lution for form factors and it does not use such con- 
cepts as "good" and "bad" current components. 

It is worth to notice that all known approaches (in- 
cluding the perturbative quantum field theory (QFT)) 
encounter difficulties while constructing a composite- 
system current operator satisfying Lorentz-covariance 
and conservation conditions ]9|-p^ . 



^ It is known that correct impulse approximation (lA) real- 
ization in the frame of traditional version of IF dynamics en- 
counters difficulties: the standard lA depends on the choice of 
the coordinate frame. We show below that lA can be formu- 
lated in an invariant way, the composite system form factors 
being defined by the one-particle currents alone. 



Similar difficulties arise in the frame of RHD ap- 
proach which is widely used in the theory of elec- 
troweak properties of composite quark and nucleon sys- 
tems §], [|l3|, At present time the FF 
dynamics is the most developed and used for composite 
systems |lo| , [r^ - [l5| , |l7| , p^ . However there are some diffi- 
culties in the FF RHD approach when the electroweak 
properties of composite systems are considered. In par- 
ticular, it was shown that the calculated elec- 
tromagnetic form factors for the systems with the total 
angular momentum J = 1 (the deuteron, the p - me- 
son) vary significantly with the rotation of the coordinate 
frame. This ambiguity is caused by the breaking of the 
so called angle condition |14 2^, that is really by the 
breaking of the rotation invariance of the theory. Some 
of the difficulties of FF dynamics are discussed in . 
A possible way to solve the problem by adding some new 
(nonphysical) form factors to the electromagnetic current 
was proposed (see and references therein). 

A different approach to the problem was proposed re- 
cently in Ref. [l3| , where a new method of construction 
of electromagnetic current operators in the frame of FF 
dynamics was given. The method of |l|l gives unam- 
biguous deuteron form factors. However, as the authors 
of note themselves, their current operator and the 
one used in Ref. |jl^ are different, since both of them 
are obtained from the free one, but in different reference 
frames, related by an interaction dependent rotation. 

Let us consider now the impulse approximation (lA) 
which is widely used for the description of composite 
systems. In I A a test particle interacts mainly with 
each component separately, that is the electromagnetic 
current of the composite system can be described in 
terms of one-particle currents. In fact, the composite- 
system current is approximated by the corresponding 
free-system current. This means that exchange currents 
are neglected, or, in other words, that there is no three- 
particle forces in the interaction of a test particle with 
constituents. It is well known that the traditional lA 
breaks the Lorentz-covariance of the composite-system 
current and the conservation law for the electromagnetic 
current (see, e.g., Q for details). 

To satisfy the conservation law in the frame of Bethe- 
Salpeter equation and quasipotential equations, for ex- 
ample, it is necessary to go beyond lA: one has to add the 
so called two-particle currents to the current operator. In 
the case of nucleon composite systems these currents are 
interpreted as meson exchange currents In the case 
of deuteron this means the simultaneous interaction of 
virtual 7- quanta with proton and neutron. However, in 
Ref. [^l[| it is shown that the current conservation law can 
be satisfied without such processes, although they con- 
tribute to the deuteron form factor. It seems that at the 
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present time there is an intention to formulate lA with 
transformed conservation properties without dynamical 
contribution of exchange currents . 

In the framework of the point form dynamics the cur- 
rent operator was constructed in the Ref. Q . The current 
operator in |^ is Lorentz-covariant and the conservation 
law is fulfilled. The approach is based on the realization 
of the Wigner-Eckart theorem for the Poincare group. 
The main idea is to extract from the current matrix ele- 
ment the relativistic invariant part - the reduced matrix 
element, i.e. the form factor, and to separate the covari- 
ant part. The form factors contain all the dynamical in- 
formation and the covariant part describes the relativistic 
transformation properties of the matrix element. 

Our approach is a generalization of the method Q for 
the case of the instant form dynamics. However, the sce- 
nario of the generalization of the Wigner-Eckart theorem 
is quite different. 

The IF of relativistic dynamics, although not widely 
used, has some advantages. The calculations can be per- 
formed in a natural straightforward way without special 
coordinates. IF is particularly convenient to discuss the 
nonrelativistic limit of relativistic results. This approach 
is obviously rotational invariant, so IF is the most suit- 
able for spin problems. 

We describe the dynamics of composite systems (the 
constituent interaction) in the frame of general RHD ax- 
iomatics. However, our approach differs from the tra- 
ditional RHD by the way of constructing of matrix ele- 
ments of local operators. In particular, our method of 
description of the electromagnetic structure of compos- 
ite systems permits the construction of current matrix 
elements satisfying the Lorentz-covariance condition and 
the current conservation law. 

To construct the current operator in the frame of IF 
RHD we use the general method of relativistic invariant 
parameterization of matrix elements of local operators 
proposed as long ago as in 1963 by Cheshkov and Shi- 
rokov ||. 

The method of ||3^ gives matrix elements of the op- 
erators of arbitrary tensor dimension (Lorentz-scalar, 
Lorentz- vector, Lorentz-tensor) in terms of a finite num- 
ber of relativistic invariant functions - form factors. The 
form factors contain all the dynamical information on the 
transitions defined by the operator. 

In the review Q| two possible variants of such kind of 
representation of matrix elements in terms of form fac- 
tors are presented - the elementary-particle parameteri- 
zation and the multipole parameterization. The variant 
of parameterization given in is an alternative one. 
In ||3^ ] the authors propose the construction of matrix 
elements in canonical basis so it can be called canonical 
parameterization. This method was developed for the 



case of composite systems in [Q|4j. The composite- 
system form factors in this approach are in general case 
the distributions (generalized functions), they are defined 
by continuous linear functionals on a space of test func- 
tions. Thus, for example, the current matrix elements 
for composite systems are functionals, generated by some 
Lorentz-covariant distributions, and the form factors are 
functionals generated by regular Lorentz-invariant gen- 
eralized functions. We demonstrate these facts below, in 
Sec. HI, using a simple model as an example. 

It is worth noting that the statement that the form 
factors of a composite system are generalized functions 
is not something exotic. This fact takes place also in the 
standard nonrelativistic potential theory (see Sec. 111(E)). 

Our formalism also gives, in fact, the description of the 
covariance properties of the operators in terms of many- 
particle as well as one-particle currents. However, the 
important feature of our formalism is the fact that form 
factors or reduced matrix elements describing the dynam- 
ics of transitions contain in the lA only the contributions 
of one-particle currents. 

So, our approach to the construction of the current 
operator includes the following main points: 

1. We extract from the current matrix element of com- 
posite system the reduced matrix elements (form factors) 
containing the dynamical information on the process. In 
general these form factors are generalized functions. 

2. Along with form factors we extract from the matrix 
element a part which defines the symmetry properties of 
the current: the transformation properties under Lorentz 
transformation, discrete symmetries, conservation laws 
etc. 

3. The physical approximations which are used to calcu- 
late the current are formulated not in terms of operators 
but in terms of form factors. 

In this paper we present the main points of our ap- 
proach. To make it transparent we consider here only 
simple systems with zero total angular momenta, so that 
technical details do not mask the essence of the method. 
We demonstrate the effectiveness of the approach by cal- 
culating the pion electroweak properties. In this case 
the canonical parameterization is very simple and can be 
realized without difficulties. The case of more compli- 
cated systems requires rather sophisticated mathematics 
for canonical parameterization of local operator matrix 
elements and will be considered elsewhere. 

The paper is organized as follows. In Sect. II we re- 
mind briefly the basic statements of RHD, especially of 
IF RHD. The IF wave functions of composite systems 
arc defined. In Sect. HI our approach to relativistic 
theory of two-particle composite systems and their elec- 
troweak properties is presented. A simple model is con- 
sidered in details: two spinless particles in the S'-state 
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of relative motion, one of the particles being uncharged. 
The electromagnetic form factor of the system is derived. 
The standard conditions for the current operator are dis- 
cussed. The modified impulse approximation (MIA) is 
proposed. The results of lA and MIA are compared. 
The nonrelativistic limit is considered. In Sect. IV the 
developed formalism is used in the case of the system of 
two particles with spins 1/2. The pion electromagnetic 
form factor and the lepton decay constant are derived. 
The model parameters are discussed and the comparison 
of the results with the experimental data is given. The 
results of calculations in lA and MIA are compared and 
are shown to differ significantly. In Sect.V the conclusion 
is given. 



II. RELATIVISTIC HAMILTONIAN DYNAMICS 

In this Section some basic equations of RHD are briefly 
reviewed. 

We use the so called instant form dynamics (IF). In 
this form the kinematic subgroup contains the generators 
of the group of rotations and translations in the three- 
dimensional Euclidean space (interaction independing 
generators) : 



J, P 



(1) 



The remaining generators are Hamiltonians (interaction 
depending) : 



N 



(2) 



The additive including of interaction into the mass 
square operator (Bakamjian-Thomas procedure ]35[ |, see, 
e.g., [Q for details) presents one of the possible technical 
ways to include interaction in the algebra of the Poincare 
group: 



U 



(3) 



Here Mq is the operator of invariant mass for the free 
system and Mj - for the system with interaction. The 
interaction operator U has to satisfy the following com- 
mutation relations: 



J, u 



VP, u 







(4) 



These constraints (^ ensure that the algebraic relations 
of Poincare group are fulfilled for interacting system. The 
relations (||) mean that the interaction potential does not 
depend on the total momentum of the system as well as 
on the projection of the total angular momentum. This 



fact is well established for a class of potentials, for ex- 
ample, for separable potentials |36|. Nevertheless, the 
conditions (^ and (Q) can be considered as the model 
ones. There exists another approach ||3^ where a poten- 
tial depends on the total momentum but that approach 
is out of scope of this paper. 

In RHD the wave function of the system of interacting 
particles is the eigenfunction of a complete set of com- 
muting operators. In IF this set is: 



Mf 



J' 



P 



(5) 



is the operator of the square of the total angular mo- 
mentum. In IF the operators , J3 , P coincide with 
those for the free system. So, in ^ only the operator 
M| depends on the interaction. 

To find the eigenfunctions for the system (|^) one has 
first to construct the adequate basis in the state space of 
composite system. In the case of two-particle system (for 
example, quark-antiquark system q q) the Hilbert space 
in RHD is the direct product of two one-particle Hilbert 
spaces: Hqq = Hq® TLq. 

As a basis in T-Cqq one can choose the following set of 
two-particle state vectors: 



\Pi ,rni; p2 ,1712) ^ \pimi) (g)\pim2) , 
{p,m\p'm') :^2pQS{p-p')Sjnm' ■ (6) 

Here pi , p2 are 3-momenta of particles, mi , TO2 — spin 
projections on the axis z, po — y'jP + M^ , M is the 
constituent mass. 

One can choose another basis where the motion of the 
two-particle center of mass is separated and where three 
operators of the set (H) are diagonal: 

\P, \/s, J, I, S, mj), 
{P, y/s, J, I, S, mj\P', Vs', J', I', S', mj') 

= NcGS^^HP-P')5iV^-^)S.IJ'SwSsS'Srn,m,, , 



NcG 



8fc Vs 



(7) 



Here = (pi + P2)(i, P^ = s, \/s is the invariant 
mass of the two-particle system, I — the orbital angu- 
lar momentum in the center-of-mass frame (C.M.S.), 
52 = (5i + 5^2)^ = S{S+l) , S — the total spin in 
C.M.S., J — the total angular momentum with the pro- 
jection TO J. 
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The basis (|^) is connected with the basis (|^) through 
the Clebsh-Gordan (CG) decomposition for the Poincare 
group (see, e.g., @): 

|P, \/s, J, /, S, nij) 



Ef dpi dp2 , ^ 
J 2piQ 2p2o 



Here 



X {pi,mi; p2,m2\P, Vs, J, I, S, mj) . (8) 



[pi ,mi; P2 ,m-2 I P, ^/s, J, S, mj] 



P2) 



The corresponding composite-particle wave function 
has the form 



(P', J', Z', m;;| 



Pc 



Nc5{P' -Pc)5jj'5nijni'j '^i'S'O^') > 



(10) 



•cO 



NcG 

4fc' 



IPc) is an eigenvector of the set (||); J(J+ 1) and mj are 
the eigenvalues of , J3 , respectively (Eqs. (^, (||)). 

The two-particle wave function of relative motion for 
equal masses and total angular momentum and total spin 
fixed is: 

iptsiHs)) = <^^Mk)k, (11) 
and the normalization condition has the form: 



J2 f ui{k)Pdk = l 



(12) 



(l/21/2mim2|5TOs)i1m,(^?,v) 



X {S I rus mi \Jmj) . 

Here A(a, b, c) = a'^+b^+c'^—2{ab+bc+ac), Yim, - a spher- 
ical harmonics, , ip are the spherical angles of the vector 
p^ (pi-p2)/2intheC.M.S., {Sms |l/2 1/2 rhi m2 ) and 
{Jmj\S Ims mi) are the CG coefhcients for the group 
SU{2), (m| D^^'^{P,p) \m) - the three-dimensional spin 
rotation matrix to be used for correct relativistic invari- 
ant spin addition. 

It is on the vectors (|^), (^ that the Poincare-group 
representation is realized in the vector state space of two 
free particles. The vector in representation is determined 
by the eigenvalues of the complete commuting set of op- 
erators: 



^P\ J\ J3 



(9) 



The parameters S and I play the role of invariant param- 
eters of degeneracy. 

As in the basis (^ the operators , J3 , P in (^ are 
diagonal, one needs to diagonalize only the operator Mj 
in (|5|) in order to obtain the system wave function. 

The eigenvalue problem for the operator Mj in the ba- 
sis (^ has the form of nonrelativistic Schrodinger equa- 
tion (see, e.g., 



Let us note that for composite quark systems one uses 
sometimes instead of the equation (n2|) the following one: 



uf{k)Pdk = 1 . 



(13) 



Here Uc - is the number of colours. The wave func- 
tion (|ll] ) coincides with that obtained by "minimal rela- 
tivization" in [Q. The normalization factors in ( [ll| ) in 
this case correspond to the relativization obtained by the 
transformation to relativistic density of states 



k^ dk 



fc2 dk 



2y/(W+W) 



(14) 



The formalism of this Section is used in the next one 
to present the method of calculation of electroweak prop- 
erties of composite systems. Particularly, the method 
of construction of electroweak current operators is de- 
scribed. 



III. THE NEW RELATIVISTIC INSTANT-FORM 
APPROACH TO THE ELECTROWEAK 
STRUCTURE OF TWO BODY COMPOSITE 
SYSTEMS 

In this Section we present our approach to elec- 
troweak properties of relativistic two-particle systems. 
To demonstrate how one describes the electromagnetic 



5 



properties of composite systems in om" version of the 
RHD instant form we first use the foUowing simple modeL 
We consider the system of two spinless particles in the S- 
state of relative motion, one particle having no charge. 
Let us note that a similar model was used in Q where the 
authors gave the description of constituent interaction in 
IF of RHD and obtained the mass spectrum. The appli- 
cation of our method in general case follows the scheme 
of this Section. The case of tt - meson is investigated in 
Sec.IV and the S =1 case in 

Electromagnetic properties of the system are deter- 
mined by the current operator matrix element. This ma- 
trix element is connected with the charge form factor 
FciQ^) as follows: 



(15) 



where p'^ , pc are 4~momenta of the composite system in 
initial and final states, = —t , = {pc—p'^Y = ^ i 9^ 
is the momentum-transfer square. The form ( p^ is de- 
fined by the Lorentz covariance and by the conservation 
law only and does not depend on the model for the in- 
ternal structure of the system. 



The Eq. (15) presents the simplest example of the ex- 
traction of a reduced matrix element, that is the sim- 
plest realization of the Wigner-Eckart theorem on the 
Poincare group. The 4-vector {pc +Pc)m describes sym- 
metry and transformation properties of the matrix el- 
ement. The reduced matrix element (the form factor) 
contains all the dynamical information on the process 
described by the current. The representation of a matrix 
element in terms of form factors often is referred to as the 
parameterization of matrix element. The scattering cross 
section for elastic scattering of electrons by a composite 
system can be expressed in terms of charge form factor 
Fc{Q^). So, form factor can be obtained from experi- 
ment and it is interesting to calculate it in a theoretical 
approach. 

In this Section we calculate the form factor of our sim- 
ple composite system using the version of RHD IF based 
on the approach of the Section II. 

Now let us list the conditions for the operator of the 
conserved electromagnetic current to be fulfilled in rela- 
tivistic case (see, e.g., [p^). 
(i) .Lorentz-covariance: 



U-\K)f{x)U{K)^Kf{K-^x). 



(16) 



Here A is the Lorentz-transformation matrix, U (A) - 
the operator of the unitary representation of the Lorentz 
group. 

(ii) . Invariance under translation: 



Here U{a) is the operator of the unitary representation 

of the translation group. 

{in) . Current conservation law: 



(18) 



In terms of matrix elements ( J^(0) ) the conservation law 
can be written in the form: 



(19) 



Here is 4-vector of the momentum transfer. 

{rv). Current- operator transformations under space-time 

reflections are: 

Up[?{x\x)J{x\x))ijp^ 

= (j\x\~x)~]{x\~x)) , 

URf{x)Up'^f{~x). (20) 

In ( pO| ) Up IS the unitary operator for the representation 
of space refiections and Up is the antiunitary operator of 
the representation of space-time reflections PT. 
{y). Cluster separability condition: If the interaction is 
switched off then the current operator becomes equal to 
the sum of the operators of one-particle currents. 
(yi).The charge is not renormalized by the interaction 
including: The electric charge of the system with in- 
teraction is equal to the sum of the constituent electric 
charges. 

In this paper the explicit equations for the form factors 
are obtained taking into account all the listed conditions. 



A. Electromagnetic properties of tlie system of free 
particles 

Let us consider first the simple two-particle system de- 
scribed in the beginning of Section HI. Electromagnetic 
current jj^^ (0) of the two-particle free system can be cal- 
culated in the representation given by the basis (^ or in 
the representation given by the basis (|^). In the first case 

the operator has the form jj^ ^ — ji^ ® I2 . Here jifj_ is the 
electromagnetic current of the charged particle and I2 is 
the unity operator in the Hilbert space of states of the 
uncharged particle. 



{pi;p2\jI°Ho)\p'i;p2] 



U~\a)f{x)ij{a)^f{x~a) 



(17) 



{P2\p'2){Pi\3nm\p'i) 



(21) 
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The matrix element of the one spinless particle current in 
the free case contains only one form factor - the charge 
form factor of the charged particle fi{Q^): 



(pi|iiM(o)|p"i> = (pi+p'i)p/i(Q' 



(22) 



So, the electromagnetic properties of the system of two 
free particles are defined by the form factor fi{Q^), 
containing all the dynamical information on elastic pro- 
cesses described by the matrix element ( |2l| ) Q . Particu- 
larly, the charge of the system is defined by the value of 
this form factor at — > 0: 



hm^ MQ^ 



/i(0) 



(23) 



Cc is the system charge. 

Now let us write the electromagnetic-current matrix 
element for the two-particle free system in the basis 
where the center-of-mass motion is separated (|^): 



(24) 



(P,V~s,\jj°HO)\P',V?) 



Here the variables which take zero values are omitted: 
J ~ S = I = 0. One can consider the matrix element 
( p^ ) as a matrix element of an irreducible tensor opera- 
tor on the Poincare group and one can use the Wigner- 
Eckart theorem, i.e. the canonical parametrization 
^3| , ||3^ giving a technical realization of this theorem. 
Thus, one can write ( p^ ) in the form 



(25) 



It is easy to understand the motivation for the parame- 
terization (|25| ) for our simple system. The 4-vector 
describes the transformation properties of the matrix ele- 
ment and the invariant function go{s, Q^, s') contains the 
dynamical information on the process. We will refer to 
go{s, Q^, s') as to free two-particle form factor. For more 
complicated systems the parameterization corresponding 
to the Wigner-Eckart theorem for the Poincare group can 
be performed using a special mathematical techniques as 
described in the papers p^, |34|, p9|]. 



So A^{s, Q , s') is defined by the current transforma- 
tion properties (the Lorentz-covariance and the conser- 
vation law): 



.s + Q')P']. (26) 



Thus, in the basis (|^) the electromagnetic properties of 
the free two-particle system are defined by the free two- 
particle form factor go{s, Q^, s'). 

So, in both representations (defined by the basis (||) as 
well as by the basis (|^)) we pass from the description of 
the system in terms of matrix elements to that in terms 
of Lorentz-invariant form factors. 

One can see that (^l|) and ( p5| ) describe electromag- 
netic properties in terms of only one form factor. Both 
of these descriptions are, certainly, equivalent from the 
physical point of view. Let us consider the difference be- 
tween these descriptions. As we will show below by direct 
calculation the free two-particle form factor go{s, Q^, s') 
is not an ordinary function but has to be considered in 
the sense of distributions in variables s , s', generated by 
a locally integrable function. So, go{s, Q^, s') is a regular 
generalizad function. All the properties of go(s,Q^,s') 
have to be considered as the properties of a functional 
given by the integral over the variables s , s' of the func- 
tion gol'S, Q^, s') multiplied by a test function. As test 
functions it is sufficient to take a large class of smooth 
functions that give the uniconvergence of the integral. In 
particular, the limit (|2^ ) giving the total charge of the 
system through two-particle form factor is now the weak 
limit: 



lim (5o(s,Q^s'),(/'(s,s')). 



(27) 



Here 0(s, s') is a function from the space of test func- 
tions. The precise definition of the functional will be 
given below. 

At the first glance it seems that the description of 
the two-particle free system in terms of the form factor 
go{s,Q'^,s') is too complicated. However, so is the real- 
ity, as we will see later in the Subsection III(E). In fact, 
this kind of description is used implicitly for a long time 
in nonrelativistic theory of composite systems, without 
calling things by their proper names. It is this kind of 
description that makes it possible to construct the elec- 
tromagnetic current operator with correct transforma- 
tion properties for interacting systems. 

The locally integrable function gois, Q^,s') can be eas- 
ily obtained by use of CG decomposition for the 
Poincare group. Using (0) we obtain for (Eq): 



(p,V^,|jW(o)|p',V7) 



dpi dp2 dpi ' dp2 ' 
2pio 2p20 2pio 2P20 



'P,Vs, I Pi;p2) 



X {pi ; P2 bl") (0) \p[ ; p1) ( p[ ; p'^ \ P' , 



(28) 
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To calculate the free two-particle form factor one has to 
use (H), (H) and the exphcit form of CG coeffi- 

cients (fel) for quantum numbers of the system. As the 
particles of the system under consideration are spinless, 
now (H) does not contain D - functions. 

It is convenient to integrate in ( p8|) using the coordi- 
nate frame with P' = , P = (0,0, P). As the result 
we obtain the following relativistic invariant form for the 
function go{s, Q^, s'): 



5o(s,Q^s') 



(S + S' + Q2)2 



(29) 



Here i?(s, Q^, s') = e{s' - si) - e{s' - S2), and 6 is the 
step function. The result, naturally, does not depend on 
the choice of the coordinate frame. 



3,2 = + 



2M- 



■(2M2 + Q2)(5_2M2) 



2M2 



The functions si,2(s, Q^) give the kinematically available 
region in the plane (s, s') (see @). 

One can see that the free two-particle form factor 
go{s,Q^,s') (29) has in fact to be interpreted in terms 
of the distributions: The ordinary limit as — > is 
zero because of the cutting d - functions and the static 
limit exists only as the weak limit (|27|). 

Let us calculate this limit. Let us define the functional 
giving regular generalized function as a functional in 
as follows: 



{go{s,Q\s'),cj>{s,s')) 
diJ.{s, s') go{s, Q^,s') (j){s, s') 



(30) 



Here 



dn{s,s') = 16 \/7^e{s- 4 M'^)e{s' -4 M'^)dfj,{s)dfj,{s') , 



(i/i(s) — -k d\/s 



(31) 



The 9 - functions in these formula give the physical re- 
gion of possible variations of the invariant mass squares 
in the initial and final states explicitly. The measure 
( ^ ) is due to the relativistic density of states (pi]), (14). 



(j){s , s') is a function from the test function space. So, for 
example, the limit of gQ{s , j s') as — > (the static 
limit) has the meaning only as the weak limit (compare 
with (El): 



lim (<?o,0) = {eA^^{s')-^l{s)),(j>) . (32) 

It is this weak limit that gives the electric charge of the 
free two-particle system. If the test functions are nor- 
malized with the relativistic density of states, then the 
r.h.s. of the Eq. ( p2| ) is equal to the total charge of the 
system. 



B. Electromagnetic structure of the system of two 
interacting particles. 

Now let us consider the electromagnetic structure of 
our simple model (p5|) in the case of interacting particles. 

As we have mentioned in Sec. II when constructing the 
bases (||) and (0) in the frame of RHD the state vector 
\pc) belongs to the direct product of two one-particle 
spaces. We can write the decomposition of this vector 
with J — I — S — mj —0 in the basis (|^). Now ( [l^ ) has 
the form: 



dP dP ' 



NcgN-cg 



x{P',V?\p',)^{Pc+p',),F,iQ'). 



(33) 



Here (P ' , vs^IPc) the wave function in the sense of the 
instant form of RHD (|l^). 
Using ( |l0|) we obtain for (^3|): 



Nc N' /— 

jj — d^/s dV s' (p{s) (p{s') 



X {Pc , V^|jp(0)|Pc ',V?)={pc+ p',)^ Pc(Q') . (34) 

We have omitted in the wave function ( pT| ) the variables 
with zero values: J = S = I —Q. 

Let us discuss the possibility of using the Wigner- 
Eckart theorem (or the canonical parametrization) in 
the case of the matrix element {pc , \/s|jp(0)|pc ' , a/s') in 
(Q). In the previous cases the state vectors and the op- 
erators entering matrix elements transformed following 
one and the same representation of the nonuniform group 
SL{2,C) 0]. Let us perform the Lorentz transformation 
of the current operator: 



(35) 
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We obtain: 



(p|?(0)b') - {p\U-\Alf{0)U{A)\p') 
= {Ap\r{0)\Ap') . 



(36) 



This means that the transformation properties of the cur- 
rent 4-vector (16) can be described using 4-momenta of 
the initial and final states, i.e. one can use the canonical 
parametrization. 

In the matrix element in the integrand of (|3^ ) the state 
vectors and the operator transform following the differ- 
ent representations of the group SL{2,C). The current 
operator describes the transitions in the system of two 
interacting particles and transforms following the repre- 
sentation with the generators of Lorentz boosts depend- 
ing on the interaction (^. The state vectors belong to 
the basis (^ and physically describe the system of two 
free particles and, so, transform following a represen- 
tation with generators which do not depend on the in- 
teraction (|^). So, if one considers the matrix element 



{pcT\/s\j^{0)\pc' TVs') (that is the interaction current 
between free states) per se , not in the context of the 
decomposition (Q), one can not use the Wigner-Eckart 
theorem. 

However, one must consider this matrix element as a 
generalized function, that is as having meaning only as 
the integrand in (^^). Let us show that in this case one 
can use the Wigner-Eckart theorem. 

The set of the free two-particle states \P , ^/s) is com- 
plete: 



dP 



NcG 

Using (|l^, (|^ we obtain 

NrNl 



dyG\p,^){p,^\ 



(37) 



d^/s dvs^ (p{k)(p{k') 



NcgN'cg 



X {p, , ^s\U-\A)j^{Q)U{A)\p, ' , Vs' 



= {pc\U-\A)iUO)IUiA)\p/) 



{Ap,\ij^{0)i\Ap,'} 



NcgNcg 
X (Ap,, VJ|j^(0)|Ap,',VI 



So, we have an analog of (|3^) in the sense of distri- 
butions and we can use the Wigner-Eckart theorem in 
the intergrand. One can speak about the Wigner-Eckart 
theorem in weak sense. 

Now the problem of canonical parametrization of the 
matrix element (^ ) can be solved if one consider the 
equality (34) as the equality of two functionals. 

Using (11), ( pl| ) we can rewrite ( p^ in the form of the 
functiona in R?T 



d^(s , s') u[k{s)) Jf^ipc ,\/s;Pc' , vs') u{k{s')) 

= {Pc+p'c)t.F,{Q^) , (39) 



J^i.{pc,Vs;pc' TVs') 



Nc N' I— 

f;—{Pc:Vs\jl,\PcWs') 



NcgN, 



CG 



The l.h.s. in ( p9[ ) contains a functional in gener- 
ated by the Lorentz-covariant function (current matrix 
element). Let us denote 



%k{s, s')^u{k{s))u{k'{s')) 



(40) 



The functional in the l.h.s. of ( p9| ) is given on the set of 
test functions ^(s , s') through an integral in and de- 
fines a Lorentz-covariant (regular) generalized function 
with the values in the Minkowski space (see, e.g., po|). 
Here is a parameter. The test-function space can be 
(in general) larger than (^). However, the uniconver- 
gence of (pa) has to be guaranteed. 

Let us write the matrix element in the form analogous 
to (H): 

J^{PctV~s;Pc\V^) = B^{s ,Q'' ,s')G{s ,Q'' ,s') . (41) 

The covariant part in (|4^) (as well as in (p5|)), the vec- 
tor Bf^{s , , s'), is supposed to be an ordinary smooth 
function and the invariant part G{s , , s') is general- 
ized function. In fact, G{s , , s') is the reduced ma- 
trix element containing the information on the process. 
This kind of representation of a Lorentz covariant gener- 
alized function as a product of a Lorentz covariant ordi- 
nary smooth function and a Lorentz invariant generalized 
function was described in [Q. 

Using ( p] ) we can rewrite (^9|) in the following form: 



d/i(s , s') V(s , s')B^{s , Q2 , s')G{s , Q2 , 5') 

= {Pc+p'c)f.F,mQ^) . 



(42) 



(38) 



To obtain the vector let us require the Eq.(p2|) to be 
covariant in the sense of distributions, that is to be valid 
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for any test function ip(s , s') in any fixed frame. The 
variation of test function in the functional ( 42 ) means in 
fact, following (^0|), the variation of the wave function of 
the internal motion. Under such a variation the vector 
in the r.h.s of ( ^ ) is unchanged as it is constructed with 
4-vectors describing the motion of the system as a whole, 
independent of the internal constituent motion. As to the 
form factor in the r.h.s it varies under the test function 
variation. So, under a variation of the test function the 
r.h.s. of ( p2[ ) remains to be coUinear to the vector {pc + 
p'c)fj.- At the same time, under arbitrary variation of the 
test function the vector in the l.h.s. in general changes 
the direction. So, for the validity of the equality ( |4^ ) 
with arbitrary test function it is sufficient to require that 
the following equation 



G{s , Q2 , s') = go{s, s') + G,{s , g2 , s') 



(45) 



(43) 



holds. This choice of the vector in (^) ensures that 
the l.h.s. of (|39| ) satisfies the condition of Lorentz covari- 
ance for the current as well as the condition of current 
conservation. 

Let us discuss the physical meaning of the representa- 
tion (|4l|), (^) for the matrix element. As this represen- 
tation is explicitly Lorentz covariant and also satisfies the 
current conservation law, then it means that the current 
operator for the composite system contains the contribu- 
tion not only of one-particle currents but of two-particle 
currents, too (see, e.g., Q): 



J 



{km) 



(44) 



Here the first term is the sum of one-particle currents and 
the second - of two-particle currents. In the case of our 
simple model each sum in (^) contains only one term. It 
is well known that if one approximates j{x) sa j'''^^ (x) 
then the current operator in IF dynamics does not satisfy 
the condition of Lorentz covariance and the conservation 
law Q . So, from the physical point of view, the covariant 
part of the current matrix element (43) which defines the 
transformation properties of the current in ( p9| ) is given 
by ( ^4| ) and contains the contributions of one- and two- 
particle currents. 

The invariant part of the decomposition (^ ) is the 
form factor or the reduced matrix element G{s , , s') 
and contains the information on the dynamics of the scat- 
tering of test particle by each of the constituents (the first 
term in (Q)), i.e. by the free two-particle system, as well 
as by two constituent simultaneously (the second term). 
So, the form factor contains the contribution of the free 
system form factor (^9|) and the contribution of some ex- 
change currents analogous to meson currents in nucleon 
systems Q. 



Here Gc is the reduced matrix element containing the 
contribution of two-particle currents (pt4|). 



Using mm, (mT), (40), mi) one can obtain from (B2 



the scalar equation of the following form: 

dy^dy^if{s)G{s,Q'^ ,s')ip{s') = FciQ^) . (46) 



The representation ( [46| ) for the charge form factor of the 
system is quite general. 

Let us note that one can use the described formalism 
in general case of composite systems with nonzero total 
angular momentum J (the detailed consideration is given 
i n ||39| ). In this case the current matrix element in the 
(|34|)-like decomposition is a matrix with matrix indices 
being the total angular momentum projections m'j , mj 
in the initial and final states. We decompose this matrix 
element in the set of the linear independent matrices : 



^^(Pc,P^)(r^(p^)p^)" , n = 0,1...2J 



(47) 



Here T^ {p'^) - is the spin 4-vector defined with the use 
of the Pauli-Lubanski vector (see, e.g., |Q, and also the 
Sec.IV(A)). 

The set of matrices ( ^ ) is the set of Lorentz scalars 
(scalars and pseudoscalars). The decomposition contains 
the 4~vectors analogous to (|4l]). 

We used the described approach to consider the sys- 
tems with J = 1 (p -meson, deuteron) and obtained good 
description of the experimental data . 

Now let us proceed with the approximate calculation 
of the form factor mt 



C. Modified impulse approximation (MIA) 

The problem of the calculation of the form factor 
G{s , , s') ( |4^ ) including exchange currents is a very 
difficult problem. We propose an approximation which 
is a kind of analog of relativistic impulse approximation. 
We propose to omit the contribution of the two-particle 
currents to the form factor G{s , , s'). 

However we will not change the covariant part of 
the current matrix element in (pl|), so that this covari- 
ant part will contain the contribution of the two-particle 
currents and so that the transformation properties of the 
matrix element will not be changed. 

So, we approximately change the generalized function 
G(s , , s') in ( ^l|) , ( ^ ) for the generalized function 
gois, Q^, s') (^, (p9|), which describes, as we have shown 
before, the electromagnetic properties of the free two- 
particle system. Nevertheless, the matrix element (IsF 
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( ^ ) as a whole will contain the contributions of two- 
particle currents, although not the full contribution but 
such that ensures its correct transformation properties. 

Let us note that our approximation does not contra- 
dict general statements (see that to obtain correct 
description of electromagnetic current of composite sys- 
tem which satisfy the Lorentz-covariance condition and 
the current conservation law one has to take into account 
many-particle currents. 

Thus, in our approximation the scalar equality (^) 
transforms into approximate scalar equality which cor- 
responds, from the physical point of view, to relativis- 
tic impulse approximation. In the developed mathemat- 
ical formalism we have not broke the Lorentz covariance 
of the current nor the current conservation law. Let us 
point out that to calculate form factor we do not use a 
special current component as it is done in other mathe- 
matical formulations of RHD (see, e.g., 10). Let us re- 
mark that, from the physical point of view, the form fac- 
tor go{s, , s') contains the contributions of one-particle 
currents only (see Equations (|2^), (p8|), (p9|)) and in this 
sense our approximation corresponds to the known im- 
pulse approximation. In order to emphasize that our ap- 
proximation differs from the usual lA we will refer to it 
as to modified impulse approximation (MIA). The form 
factor of the composite system in MIA has the form: 



Fc{Q^) = J d^sdV?^is)go{s,Q\s')ip{s') 



(48) 



We do not discuss in this paper the problem of going 
beyond the limits of MIA and of obtaining corrections to 
5o(s ,s')iii (p5|), (p^). This means that if considering, 
for example, nucleon systems we do not take into account 
meson current. 

Let us consider now the fulfilling of the conditions (i)- 
(vi) for the electromagnetic current. 

The conditions (i)-(iii) are satisfied by construction. 
For example the fulfilling of (i) and (iii) is ensured by 
the correct transformation properties of the 4-vectors in 
(H), (|ll), and (H). 

The condition (iv) is satisfied immediately as the 
form factor go{s , s') in ( p5[ ) and the form factor 
G{s , , s') in (^ ) are scalars in our simple model ^ 

The condition of cluster separability (v) needs a more 
detailed consideration. At large distances (or if the 



The currents which do not conserve the parity also can be 
considered in our formalism. In that case one can separate 
not only the scalar part of current matrix element but the 
pseudoscalar part, too. This case is considered elsewhere. 




interaction is switched off) the contribution of two- 
particle currents has to go to zero: G, 
in (p5|). This means that in the form 
tor G{s , , s') has to transform into go{s , , s'). Let 
us remark that the condition of cluster separability is 
fulfilled in MIA, too, as in this approximation the use 
of go{s , , s') instead of G{s , , s') is supposed from 
the very beginning. When the interaction is switched off 
the generalized function go{s , s') for the free two- 
particle system acts on a larger space of test functions 
than (|40|). As go(s , , s') contains only the one-particle 
current contributions (Eq) the condition (v) is satisfied 
and the composite-system current go to the sum of the 
one-particle currents. 

The condition on the charge to be nonrenormalizable 
also is fulfilled directly in MIA because the weak limit 
( p^ ) does exist on test functions (|40|). 

So, our prescription for the construction of the cur- 
rent in MIA satisfies all the conditions for the current 
operator. 

Let us note that the equation ( ^8| ) for the composite- 
system form factor is analogous to the equations obtained 
in the framework of the dispersion approach ^3 41 1 (see 
also [Q) based on the analytic properties of the scat- 
tering amplitudes, matrix elements, form factors in the 
complex energy plane. 

As the dispersion approach is rather correctly derived 
in the frame of QFT , this fact can be considered as a 
possible link between QFT and RHD. The establishment 
of such a link is one of the unsolved problems of RHD 

i- 

Let us note that an immediate application of the ap- 
proach to quark systems is difficult to realize because of 
the fact of quark confinement. However, there are some 
investigations based on similar ideas where the form fac- 
tors of hadrons as constituent -quark bound states are 
considered in the frame of the dispersion technique of 
the integral over composite particle mass Es] . 



D. MIA versus I A 

Let us compare the approximation MIA with the well 
known I A. 

To do this let us first calculate the form factor in IF 
RHD not using the canonical parameterization. In par- 
ticular, let us formulate I A in terms of operators as it is 
formulated usually (not in terms of form factors). Let us 
decompose the matrix element ( p^ ) through the complete 
set of states (O): 



[Pc\UO)\Pc 



dpi dp2 dpi ' dp2 ' 
2pio 2p20 2p' 2p' 
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x(Pc \pi;p2){pi;P2 \jtJ.\pi';P2') 
X {pi';p2'\p'c) ■ 



(49) 



Here {pi]p2 \pc) is wave function of constituents in com- 
posite system. 

If the current matrix element in (^9|) is taken in the lA 
approximation (Q) and contains one-particle currents 
only, then the Eq. ( ^ ) is sclfcontradicting Q . 

To write the form factor in terms of wave functions 
( p^ ) one has to perform the CG decomposition of the 
basis in terms of the basis in the wave functions 
( ^ ) and to use the exphcit form for CG coefficients (||) 
for the quantum numbers of the system: 



[Pl]P2 \Pc) = ^J^{P ,y/s ,J ,1 ,S ,mj \pc 



(50) 



The current matrix element in (|49|) has the form ( 21 ) . 
The one-particle currents are expressed through the form 
factors 

The Eq.(^9|) is an equality for two 4- vectors. Taking 
different components of this equaiity and expioiting S- 
functions in integrais, one can caicuiate the form factor 
of the composite system. The resuit of caicuiation of the 
form factor in this way is not unambiguous. In particuiar, 
it depends on the actual choice of the component of the 
current ( |49| ) to be used in the caicuiation. Moreover, the 
resuit depends on the coordinate frame chosen to perform 
the integration in (|4^) . This is the general feature of lA 
in the usuai formulation of IF RHD (see, e.g., [Q). 

Let us write the final result of the caicuiation of the 
form factor from the equation for the nuil-component of 
the current and performing the integration in the coordi- 
nate frame where Pc ' = , Pc = (0, 0,p). If now we write 
the integral in terms of the invariant variabies s, s' the 
obtained form factor has the form: 



Mc V2(2M2 + g2) 



4Af2-f Q2 



dy/s dys 



(s + s' + QYQ^ 



s' ^(s -4M2)(s' -4M2) [A( 
1 0(8, Q\s') 



(ss')i/4 V^(s + Q2) 



Q2,S')]'/' 

^(s)^(s')/i(Q')- (51) 



The Eq.(|l|) differs from (H), obtained with the use of 
the two-particle free form factor. In the case of wave 
fmictions satisfying the conditions (|ll|), (p^), the form 
factor ( ^ ) satisfies the normaiization: Fc{0) = Cc- Let 
us note that the form factor obtained in this way from 



the third current component in ( |49[ ) docs not satisfy this 
condition. 

Let us compare lA and MIA results. Let us note once 
again that in MIA we separate (by use of the scheme of 
canonical parameterization) the covariant part of the cur- 
rent matrix element in (|4^) prior to perform any calcu- 
lations. This covariant part ensures the correct transfor- 
mation properties of the corresponding decompositions 
in terms of free-particle states. The difference between 
(H) and (|5l|) is: 

AFe(g2)= J dy/ldV? f{s)^{s') 



xgo{s,Q^,s') [l-i^(s,Q^s')] 



(52) 



Ris,Q^,s') = 



Mc ^2 (2 j\f 2 + Q2 
2 4 Af2 + Q2 



(S + S' + Q2)S 



1 



(53) 



The value R{s, Q^, s') presents an additional factor to 
one-particle currents, that is in reality the two-particle 
current contributions. This term ensures the Lorentz co- 
variance of the electromagnetic current matrix element 
and the current conservation law in ( |39| ) . Let us note that 
this additional term contains no dynamical information 
on the interaction of test particle with two constituents 
simultaneously. It does not depend, for example, on the 
interaction constants for such a process. 

So, to summarize, we can write the following schematic 
equations: 

{IA)Brert + {IA)Lab 



(MIA) 



Breit 



{MlA)Lab 



It is well known that the standard lA depends strongly on 
the coordinate frame used for the calculation. The MIA 
results do not depend on it at all. So, the differences 
between lA and MIA results for different lA coordinate 
frame can be rather significant. 

Notice that lA and MIA coincide in the nonrelativis- 
tic limit. As this takes place, the nonrelativistic limits 
of form factors, which were obtained from the different 
current components, are identical. Hence the difference 
between lA and MIA is connected with the breaking of 
relativistic covariance conditions really. 

We give the quantitative comparison of the form fac- 
tors obtained in lA and in MIA in the Section IV where 
the realistic calculation of the pion electromagnetic struc- 
ture is given. 
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E. The nonrelativistic limit 

The description of composite-system form factors in 
terms of distributions is not a specific feature of our rel- 
ativistic approach. The similar formalism is used in non- 
relativistic theory of composite systems for a rather 
long time (although not referring to the mathematics of 
distributions). In the nonrelativistic limit our approach 
gives the formalism developed in ]47| . 

In the nonrelativistic limit the relativistic charge form 
factor MS) has the following form: 



Fnr{Q^)^ / kUkk'^dk'u{k)goNR{k,Q^,k')u{k') 



goNR{k,Q'^,k') 



e{k,Q^,k') = dik'- 



kk'Q 



k-^ 
2 



e{k, 



,k') 



(54) 



(55) 



-dik' 



" 2 



Here goNR^I^jQ'^ , k') is the free relativistic form factor 
obtained from ( p9| ) in the nonrelativistic limit. /i(Q^) 
is the charged-particle form factor. The obtained result 
coincides with that derived in standard nonrelativistic 
calculations |^ . 

Rigorously speaking, the Eq. ( |5^ ) has to be interpreted 
as a functional in the sense of distributions generated by 
the function goNR{k, Q^, k') and defined on test functions 
u{k)u{k'). The ordinary function ( [55[ ) generates regular 
generalized function defined generally on the larger class 
of test functions ip{k,k') in R^, providing the uniform 
convergence of the integral. One needs the uniform con- 
vergence to take limits in the integrands. 

Let us define the functional in by the following 
regular distribution (compare with ([3C|)-(pl|)): 



(50JVi^(fc,Q^fc'), V'(fc,fc')> 
dn{k, k') goNRik, Q'^,k') ip{k, k') , 



(56) 



d/i(fc, k') = e{k) e{k')dfj,{k) d^i{k') , dfi{k) ^k^dk . 

The function goNRik,Q^ , k') which appears in |47| ] 
quite formally, here has a definite physical meaning 
and describes the electromagnetic properties of non- 
relativistic free system of two spinless particles in the 
S - state, one of particle having no charge (compare 
with go(s,Q^s') in (||), (|9|), (^)). The statical limit 



limQ2_,Q gaNR{k,Q'^ ,k') giving the system charge exists 
only in the weak sense as the limit of the functional (|5^): 

lim {gaNR{k,Q'^,k') , %l;(k,k')) 



- {e,5{^i{k')-^l{k)),^^{k,k')) 



(57) 



On the test functions ?/;(fc, k') = u{k) u{k') (with u{k) 
- being the normalized bound state wave function), the 
functional ( |56| ) defines the bound state form factor in the 
nonrelativistic lA (jsj). The weak limit (|5^) is equal to 
the system charge: 



lim (goNRik,Q'^,k') , i}{k,k')) 



k^ dk v?'{k) ~ Ec 



(58) 



To go beyond nonrelativistic lA one has to addend 
some terms to goNR{k, Q^, k'). For example, such terms 
cause the meson exchange currents in two-nucleon sys- 
tems. So, in the standard nonrelativistic theory the dy- 
namical treatment of exchange currents is performed in 
the same way as in our relativistic approach (psf). 

So, to conclude, one can consider our approach to lA to 
be a relativistic generalization of nonrelativistic lA, and 
our equations for form factors in this approximation to be 
a relativistic generalization of the equations of |^ . Let 
us remark that in more complicated systems (e.g., in p 
- meson and deuteron) our relativistic form factors also 
have correct nonrelativistic limits which coincide with 



IV. THE ELECTROWEAK STRUCTURE OF 
PION 

Now we apply the method of previous sections to the 
calculation of the electroweak structure of pion. There 
exists a lot of experimental data on pion, so the effec- 
tiveness of the method can be checked by the comparison 
with the data (see, e.g., |13] and references therein). 



A. The electromagnetic form factor of pion 

The pion is spinless, so the electromagnetic current ma- 
trix element has the form ( p^ ) with Pc ^ Ptt , Fc{Q^) 
Ptt{Q'^)- In the frame of composite quark model pion is 
considered as the bound state of u- and d- quarks. We 
suppose that quark masses are equal: m„ = = M . 

To calculate in MIA the composite-system form factor 
one needs to construct first the free two-particle form 
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factor (p5[), (p9|), ([48|). Contrary to the simple model of 
the previous Section now we consider the system of two 
charged particles with spins 1/2. This gives the follow- 
ing complications. First, the Eq. ( pT|) for the current 
operator of the free system is now transformed to the 
form: 



(59) 



Here J(i,2)/j - the electromagnetic currents of particles, 
^(1,2) ~ the unity operators in the one-particle state 



Hilbert spaces. The Eq.(59) can be rewritten in terms 
of matrix elements: 

(pi , mi ; p2 , TO2 14"^ (0) \p[ , m[ ; p!, , ma) = 

= {P2,m2\p'2,m'^){pi,mi\jif,\p[,m[) + (1^2). (60) 

Second, the matrix element of one-particle current 
contains now, contrary to (p2|), the magnetic form factors 
of quarks as well as the charge ones. Now the parame- 
terization (the elementary-particle one following jij) is 
of the form: 

{p, m\j''{0)\p', m) 

= Upml'^Up'.n' Fi{Q^) - UprnCr^''qu Wp'm' ^2(Q^) , 

(61) 

Uprn - the Dirac bispinor, - Dirac matrix, 

t'-" = lirY -l^'in , q.^ip-p%, 

Using multipole parameterization we can write the one- 
particle current matrix element in terms of Sachs form 
factors: 



4M^ 



Gm{Q^)=F,{Q^) + kF2{Q^) , 

Fi(g^)=eA(g^), F2{t) ^ ^F2{Q') . (62) 

Here Ge.m - Sachs electric and magnetic form factors, 
respectively, e is the particle charge, k, is the anomalous 
magnetic moment. 

It is convenient to use the canonical parameterization 
of matrix elements : 

(P, j^(0) m') 



E 



m\D={p, p')\m"){m"\h{Q^)K+^f2{Q^)R^m') 



^; = b + ?''V, Rt^^e^uXflP^p'^-TP^p') . (63) 
r(p) is 4-vector of spin: 

Am) 



T{p)^M] + 



Po + M 



ro(p) {pj) ■ 



The form factors fiiQ^), f2{Q^) are the electric and mag- 
netic form factors of particles. They are connected with 
Sachs form factors UM: 



2M 



4 



■■Ge{Q') , 



M^4Af2 



GM(g' 



(64) 



Third, now the CG coefficients are of more complicated 
form. They are given by (^ with J — S — I = 0. Con- 
trary to the previous simple case, now the CG coefficients 
contain the Wigner rotation matrices. 

Finally, the free two-particle form factor for the system 
of two particles with spin 1/2 and quantum numbers J — 
5 = Z = is of the form (see also [|9|): 



So'(s,Q^,s') = ric 



{s + s' + Q^)Q^ 



2^(s-4M2)(s'-4Af2) 



[A(5, -Q2, s')]3/2 ^1 + Q2/4^p 

X {(s + s' + g2)(G^(g2) + G|(g2)) cos(wi + lo2)+ 



+ j^as, s')(Gl,{Q^) + GUQ^)) sin(c.i + c.^) 



(65) 



Here 



e(s, g2, s') = v/ss'g2-M2A(s,-g2,s') , 

ric is the number of quark colours, loi uj2 - the Wigner 
rotation parameters: 



uJi = arctan ■ 



C(s,g2,s'^ 



M 



(jj2 — arctan 



(V^ + vV)2 + g2j -I- Vss'(V^ + V7) 
a{s,s')i{s,Q^s') 



M{s + s' + g2)a(s, s') + Vss'(4Af2 4, q2) 

(66) 
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with a(s,s') =2M + ^+^/7, and G^j^jiQ^) 

are Sachs 

form factors for quarks. The - function in (|6^) is the 
same as in (p9|). 

An interesting effect follows from (^5|): due to the rel- 
ativistic Wigner spin rotation effect the pion charge form 
factor contains the contribution of quark magnetic form 
factors. 

The pion charge form factor can be calculated using 
with (|6^) for the free two-particle form factor: 



dV^dVs' <f{s) 5««(s , Q2 ^ s'Ms'). (67) 



B. The lepton decay constant of pion 

The lepton decay constant is defined by the 
electroweak-current matrix element [[15|: 



1 



(27r)3/2 



(68) 



Pti - 4-momentum of meson. Let us decompose the l.h.s. 
of ( |68| ) in the basis ( fzl). Using the explicit form of the 
meson wave function (noh one can obtain for (68): 



NCG 



(27r)3/2 



(69) 



As in Section II (Eq.([4l|)) one can divide the integrand in 
( |69| ) into two parts: the covariant part (smooth ordinary 
function) and the invariant part. 



CG 



{0\jAO)\p^,V~s)=iGis)BAs) 



1 



(27r)3/2 



(70) 



The invariant form factor G{s) is a gen eralized func- 
tion. In the same way as in calculating ( [46| ) of the previ- 
ous section, we now obtain the lepton decay constant of 
pion in the form 



dy/s G{s)ip{s) = U 



(71) 



In general, the form factor G{s) can be calculated in 
the frame of the standard model for electroweak inter- 
actions. However, in this paper we limit ourselves by 
4-fermion interaction. We take for G{s) the form factor 
which parameterizes the decay of free two-quark system: 



(0|jr(0)|P,Vi)-*Go(.)P^ 



(27r)3/2 



(72) 



The explicit form ([72|) is written by analogy to (|2^) not 
taking into account the current conservation law. The 
form (|7|) is quite similar to (^ but instead of the con- 
stant the form factor depending on invariant variables 
is written. To calculate G'o(s) let us decompose ( |7^ ) in 
the one-particle basis (||). Now we obtain for ([72|): 



E 



iGois)P^ 



dpi dp2 



1 



(27r)3/2 



2pio 2p: 
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mi ,7712 ,2c 

X {Pi ,mi;p2 ,m2\P ,y/s}- 



(73) 



ic = 1, 2, 3, the sum over ic is the sum over the colours. 
The CG coefficients are known (^. The current matrix 
element in the basis (^ can be written in the standard 
way in terms of the lepton decay current matrix element 
ilBl: 



i2n) 



(Obl°^|Pi ,mi;p2,m2) 
■ v{P2 , "^2)7^^(l + l^)u{pi , mi) . 



(74) 



We integrate in ( |73| ) in the coordinate frame with P = 0. 
Finally, we obtain: 



Go(s) = 



2V2 7rPo 



(Po + M) 



1 - 



(Po + Mf 



(75) 



Substituting ( [75D in the Eq.(|7l|) we obtain the result 
which has the following form if written in invariant vari- 
ables: 



2Mn, /• , ^ 1 , . 



(76) 



Let us notice that the Eq.([76|) coincides with that ob- 
tained in the frame of light-front dynamics |l^. How- 
ever, although all forms of RHD are unitary equivalent 
|l2| , nevertheless after the physical approximations are 
made in more complicated cases the results, e.g. for 
form factors, can be different. This is possibly due to 
the fact that the unitary operators connecting different 
forms of RHD are interaction dependent ||l^ and so the 
RHD forms realize one and the same approximation in 
different ways. 

Let us remark that the nonrelativistic limit of the 
Eq.([76|) gives the standard form in terms of coordinate 
space wave function at zero value. 
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C. The results of calculations 

To calculate the electroweak structure ol pion using 
(|67|), (65), (1?^), ( [Ill ) the following meson wave functions 
were utilized: 

1. A gaussian or harmonic oscillator (HO) wave func- 
tion 



u{k) = NHoexp{^k^/2b^) 
2. A power-law (PL) wave function 
u{k) ^ NpL{k^/b^ + n: 



2,3 



(77) 



(78) 



3. The wave function with linear confinement from 
Ref. [El: 



u{r) = Nt exp(-ar^/^ - /3r) , 



;VMa, 



M 



(79) 



a,b- parameters of linear and Coulomb parts of potential 
respectively. 

In the Ref. ||l^ in the calculation of pion electromag- 
netic structure we supposed the quarks to be point-like. 
The results of JTot can be considered as preliminary re- 
sults. However, one has to take into account the structure 
of constituent quarks po|| , in particular, the anomalous 
magnetic moment. As anomalous magnetic moments are 
connected with finite size of quark, one has to take into 
account the explicit form of quark form factors entering 
( |65| ) and the pion charge form factor (|6^). As in ||l^ let 
us use the following forms for quark form factors: 



(80) 



Here Cq - the quark charge, 



the quark anoma- 



lous magnetic moment (in natural units). To obtain 
the explicit form of the function f{Q^) let us consider 
the asymptotics of pion charge form factor as — > 
oo , M ^0. 

To obtain the asymptotic behavior let us first make 
the asymptotic estimation of the integrals in (^) in the 
point-like quark approximation (/(Q^) = 1 , k — 
indSOj) ). Omitting the details of calculation (given in 
we write the final result for the asymptotics in the 



form: 



(81) 



The asymptotics does not depend on the actual form of 
the wave function and coincides with that obtained in 
QCD. The actual form we obtain, e.g. for (77) is: 



F^iQ^) -32V2. 



It is worth to compare the form 
QCD resuh g: 



(82) 



82) with the detailed 



(83) 



lias/n ~ 0.1 then (H) and (H) coincide at & - 0.1. 
So the asymptotics ( |8l] ) is quite reahstic. 

In the case of non-point-like quarks we obtain another 
asymptotics because the form factor depends upon the 
momentum transfer. It is known that QCD gives log- 
arithmic corrections to (|8l|). To agree with this QCD 
corrected asymptotics we can, for example, choose the 
following form for f{Q^): 



im = 



1 



1 + ln(l + (r2)Q2/6) 



(84) 



Here (r^) is the MSR of the constituent quark which can 
be considered as the model parameter. Let us fix it (as 



in |18|) to be: (r^) ~ O.S/M^. 



For the constituent quark mass in pion we use the value 
which is usually used in the calculations in RHD: M = 
0.25 GeV. 

The quark anomalous magnetic moments can be taken 
from 0: k„ = 0.029 , Kd = - 0.059. 

We choose the parameters h in (^, (^ and a in 
( [79| ) in such a way as to fit the pion MSR: (r^) — 
(0.432 ± 0.016 ) Fm^ |53). We choose this way to fix 
the model parameters because the pion MSR is defined 
by the form factor at small values of , that is the range 
where potential models work well. 

The fit of the pion MSR gives the following parameters 
of the wave functions: in the model ( |77| ) h ~ 0.2784 GeV, 
model (|7|) at n = 2 6 = 0.3394 GeV, model (^) at n = 
3 & = 0.5150 GeV, model (H) 6 = (4/3)a^, Us = 0.59 at 
light meson mass scale, a = 0.0567 GeV^. 

The results of calculation are presented on Figs.l and 

2. 
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FIG. 1. The square of the pion form factor at small values 
of momentum transfers for different models. 

The square of the pion form factor at small values of 
momentum transfers for different models (|7^) - (|7|) is 
presented on Fig.l. Results of calculation in the models 
(H), (H) at n = 3 and (H) coincide very closely. 

The calculations of product Q'^ F.„{Q^) at high mo- 
mentum transfers for different models (|7^) - ( [79| ) are 
presented on Fig. 2. Legend is following: 1 - harmonic 
oscillator wave function (|7^), 2 - power-law wave func- 
tion (78) at 71 — 2, 3 - power-law wave function ( [ts] ) at 



n = 3, and wave function from model with linear con- 
finement (^9|) (these curves coincide very closely). 

All the models for the interaction (|7^), (^), (^ give 
a good description of the existing experimental data 0. 

The dependence of the results on the actual model 
is much less pronounced that in the case of point-like 
quarks |]l9| . 

The lepton decay constants calculated following Eq. 
( [tgI ) with different wave functions have the following val- 
ues: ^ 0.1210 GeV in the model = 0.1327 
GeV in the model (|7|) with ?i = 2, = 0.1282 GeV in 
the model (||) with 7i = 3, and = 0.1290 GeV in the 
model (^^. Let us emphasize that we have used no fit- 
ting parameters to calculate the lepton decay constant. 
Nevertheless, the obtained values are very close to the 



experimental value: /^r < 



^The JLab new results 
our approach in Q 



0.1317 ± 0.0002 GeV [M 



are discussed in connection with 



0.0 



Q'(GeV ') 

FIG. 2. Electromagnetic form factor, f,r(Q^), at high 
momentum transfers. 



Now let us compare the numerical results for the pion 



form factor obtained in MIA (67) with that of the tradi 



tional lA. Let us choose for the comparison, for example, 
the null-component of the current. 

To obtain the pion form factor in lA we proceed in 
the same way as while obtaining ( pT| ) of the preceding 
Section. Now, however, 

1) the decomposition ( |l5| ) of the I A matrix current 
element over the state set (^ is realized following (|60|), 

2) the parameterization of the one-particle matrix el- 
ement is given by (|6|), (^) (instead of (p2|)), 

3) the GG coefficient (^ in (|5^) are for pion quantum 
numbers. 

Acting in the same way as while obtaining (|5^), and 
using the null-component of the current matrix element, 
we can write the pion form factor in lA in the following 
form: 



M^V2(2M2 + Q2) 

As" d^/s d\/s' 

1^ VTs^4M2)(s'-4l?2) 



(p{s)(p{s') 



X {(.S + S' + Q2) [g« (Q2) + G|(Q")] COs{lOi + U2) 

+ l^^{s,Q',s') [GliiQ') + Gii{Q^)\ sin(c.i+c.2)}. 



(85) 



17 



Here = 139.568±0.001 MeV |4| is mass of pion. 

The normalization condition -F'7r(0) = 1 is satisfied for 
the form factor ( |85| ) if the wave functions ( pi] ) satisfy 

To compare the numerical results given by the 
Eqs.(|67|), ( pSi ) with that given by (^) let us calculate the 
pion form factor using the wave function ( [tt] ) with the 
parameters of the calculations presented in Figs.l and 2. 
The results are shown in the Fig. 3. The results obtained 
with the use of the parameterization (Esl), (65) differ es- 



sentially from that obtained without such parameteriza- 
tion (85). The form factor calculated in our approach 



describes the existing experimental data adequately. 




0.0 



Q'(GeV') 

FIG. 3. QV(Q2) for MIA mand for lA (2). Results of 
calculation with wave function ([77|). Parameters are the same 
as in Fig.l. 

Let us emphasize once again that the form factor ob- 
tained in MIA does not depend on the choice of coor- 
dinate frame. This is an important advantage of our 
relativistic MIA. 



magnetic current also the conservation law auto- 
matically. 

The method of the construction of the current op- 
erator matrix element consists of the extraction of 
the invariant part - the reduced matrix element on 
the Lorentz group (form factor) - and the covariant 
part defining the transformation properties of the 
current. The form factors contain all the dynamical 
information about transition. The properties of the 
system as well as the approximations used are for- 
mulated in terms of form factors, which in general 
have to be considered as generalized functions. 

The approach makes it possible to formulate rela- 
tivistic impulse approximation (modified impulse 
approximation - MIA) in such a way that the 
Lorentz-covariance of the current is ensured. In 
the electromagnetic case the current conservation 
law is ensured, too. 

The results of the calculations are unambiguous: 
they do not depend on the choice of the coordinate 
frame and on the choice of "good" components of 
the current as it takes place in the standard form 
of light-front dynamics. 

The effectiveness of the approach is demonstrated 
by the calculation of the electroweak structure of 
the pion. Our approach gives good results for the 
pion electromagnetic form factor in the whole range 
of momentum transfers available for experiments at 
present time. 
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V. CONCLUSION 



Let us summarize the results. 

1. A new approach to the electromagnetic proper- 
ties of two-particle composite systems is developed. 
The approach is based on IF RHD. 

2. The main novel feature of this approach is the new 
method of construction of the matrix element of 
the electroweak current operator. The electroweak 
current matrix element satisfies the relativistic co- 
variance conditions and in the case of the electro- 
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